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' Linear response methods applied to electron systems often display a level of accuracy which 

^\ I is notable when viewed in terms of the strengths of perturbing interactions. Neglect of higher 

response terms is in fact justifiable in many cases and it can be shown to stem from an intrinsic 
5^ ' interference between atomic and electronic length scales. For fluid metallic systems it can be further 

^1^, shown that electron-ion structure (increasingly accessible experimentally) can be understood from 

■ an application of linear response in the electron system, combined with hard-sphere like correlation 

for the ionic component. 
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^ ' The nearly-free-electron (NFE) approximation underlies much of our understanding of the properties of condensed 
C/3 ' matter, in particular simple metals. While ab-initio simulation techniques have long superseded the NFE approxi- 
mation in quantitative accuracy, it remains an important source of insight and of simplifying concepts to elucidate 
qualitative trends across different materials. It also provides guidance in situations that remain out of the reach of 
computational ab-initio techniques For many years, the density p*"'^(k) of an initially uniform electron gas induced 
by an embedded pseudo-potential v^^ (k) has been successfully treated at linear order even though v'^^ (k) is not neces- 

■ sarily a small perturbation. The linear approach is a key component in many applications of the NFE approximation, 
\ examples of which include pseudo-potential calculations of the free-energy of simple metals, their relative structural 

stability (and corresponding cohesive properties) and also the determination of effective ion-ion potentials |^|^. The 
accuracy of the latter is a particularly striking example of the efficacy of linear response; while the energy scale of 
Q unscreened ions at typical separations is of the order of Ry, linear screening leads to ion-ion potentials fully capable 
^ of describing observed structural phase transitions and implying consequent energy scales of the order of mRy. 

, Here we address the evident success of the linear approximation which to date remains incompletely resolved. We 
^ ■ show that the implied neglect of higher order response is supported by physical arguments. In particular we explicitly 
^ \ demonstrate that the nonlinear terms are small for specific cases, and give arguments to suggest that this may be 
• expected to hold more generally, the main exception being hydrogen. As an application of the underlying argument, 
CnJ ' but one with experimental consequences, we demonstrate that simple linear-response theory augmented by a hard- 
en sphere approximation for ionic structure leads to a quantitatively accurate analytical representation of electron-ion 

■ structure factors Sei{k) in liquid metals, these now in principle accessible through recent advances in both neutron 
and x-ray scattering techniques. Another route to effective electron-ion interactions therefore opens, but here through 
the fluid state. 

To begin, consider the response of the interacting electron gas to a single ion, where the electron-ion interac- 
"j^ ■ tion is modeled by a pseudo-potential, taken as a simple local one- parameter empty-core form |^ i.e.: v^'^{k) = 
—{A'Ke^/k^)cos{kRc), where the atomic core-radius Rc or equivalently the zero-crossing fco = ■k/2Rc is typically fixed 
I by an atomic property such as the ionization energy (or by a measurable crystalline metallic property such as the Fermi 
surface 0). The pseudo-potential leads to a local electron density inhomogeneity representable by p™'^(k). There 
Ch , are two routes to represent this induced density, the first (esentially exact) from solving the familiar self-consistent 
^ Kohn-Sham equations within the local density approximation (LDA) and the second from the standard expansion 
of the response in powers of the perturbing (pseudo)potential, i.e.; 

r\k) = Xi{k)vP^k) X2{k, k,, k2)vP^{k^)vP^{k2) + .... (1) 
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Here the response functions Xn{ki...) are properties of the homogeneous interacting electron gas, the first being the 
well known linear response function . The second is given by: 



X2(fci, fc2, fca) = (x°2{ki, k2,k^) + ]^iji2(k,, k2,k3)x°Aki)Xi{k2)Xi{k3?j I (e(fci)e(fc2)e(fc3)) , (2) 

where e(/c) is the usual dielectric function: e(fc) = 1 — _f47re^//j^ /ii(fc)) Xi(fc)- In (2) the x°(A:i...) are the non- 
interacting response functions (known to second order H) and the /i„(/ci..fci+i) are the homogeneous limits of the 
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n-th functional derivatives of the exchange-correlation potential with respect to density. In particular, /ii(fc) is related 
to the spin-averaged local field correction (LFC), G(fc) = (fc^/47re^)^i(fc). Fig. 1 compares the full LDA response, and 
equation (||), taken to second order; note that it appears to capture most of the complete response with considerable 
accuracy. 



Cl 

2" 




-0.020 



k/2k. 



FIG. 1. A comparison of full non-linear LDA response [p{k)—p'^^\k)] (solid line) to second order LDA response (dashed 
line) for an empty core pseudo-potential with Rc=l-5ao embedded in an electron gas with density parameter rs—Sao- For the 
scale, compare this to the full response with the limit p(fc=0) = l. The higher order response is of the order of a few % of 
the full response and in turn, the second order response captures almost all the non-linear response. (The small difference at 
fc ^ is a numerical artifact stemming from the use of a large but finite real-space cut-off radius utilized in the Kohn-Sham 
procedure.) In the insert is plotted the maximum of the 2nd order response vs. Rc/rs for rs=2ao (dotted), rs=3ao(solid) and 
rs=5ao(dashed). Note especially the minimum at Rc/ra=0.41 which corresponds to ko=2kF- It is reduced by an order of 
magnitude from the value at Rc=0 (hydrogen) and is traced to an interference between atomic and electronic length-scales. 

Interestingly enough, the combined effects of exchange and correlation partially cancel between first and second 
order, and this implies that the neglect of the /x„(fci..fci+i) at both orders (the RPA) is found to be more accurate 
than the result obtained by merely including them at a single order only [ p^ . This has important implications for the 
widespread application of linear response theory in the derivation of effective ion-ion potentials in (simple) metals; 
the neglect of higher order response results in an overestimate of the role of exchange and correlation. The accuracy 
of the second order response depicted in Fig. 1 also implies that the use of more accurate LFC's could, in some 
cases, lead to an improvement in accuracy over a full Kohn-Sham LDA calculation. Although much effort has gone 
into obtaining LFC's beyond the (fc = 0) LDA limit at linear order the second order LFC, directly related to 
/i2(fci, ^2, ^3), to date remains unknown beyond the LDA form. However, the second order electron LFC is the direct 
analog of the third order direct correlation function c*^'^-'(fci, fc2, A:3) of classical liquid-state theory for which various 
successful approximations based on lower order correlations functions have been derived [p^ . (Since the electron 
liquids are more weakly correlated than their classical counterparts p3| , it might now be suggested that application 
of these classically inspired approaches to the electronic case would be useful.) 

A central question now arises (whose answer is important to the proposition we make on electron-ion structure): 
Why is the non-linear response contribution depicted in Fig. 1 evidently so small? An immediate possibility is that 
higher order terms in equation (|^) are large, but actually vary in sign and therefore mutually cancel, order by order. 
But another is that the higher order terms are each individually very small. The success at the level of 2nd order 
response evidently implies that the latter is the case: we find that the response series converges very rapidly. This 
might be physically anticipated since a larger atomic-parameter Rc implies a smaller perturbing potential, and the 
non-linear response shows a clear decline with increasing Rc. As anticipated the second order response is found to be 
largest for Rc — (hydrogen), but as Rc increases from zero a noticable secondary minimum occurs when the inverse 
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atomic length kg is equal to 2kp. For the cases plotted in Fig. 1, the secondary minimum is reduced by an entire order 
of magnitude when compared with the value calculated for hydrogen, and is typically factor of three lower than the 
secondary maximum at larger Rc. This minimum is attributed to the following; the second order response function, 
X2{k, ki, ^2), itself peaks when the summed arguments in (|l]) are close to 2kp ||lo|. Accordingly, if the pseudo-potential 
zero-crossing fco is near the response peaks at 2kp, a maximal cancellation or maximal destructive interference of the 
atomic and electronic lengthscales occurs, leading to a minimum in second order response. We may now postulate 
that for the simple metals a similar interference effect occurs for the higher order terms of (|^). 

Typically the value of kQ/2kF lies between 0.75 and 1, and is therefore very close to the secondary minimum in 
the non-linear response. Note that the ratio of the atomic and electronic length scales is set primarily by the volume 
energy terms in the total ground state energy, and is almost independent of structure |^ . This clarifies in large part 
why the ubiquitous linear response approximation performs so remarkably well for many materials and why the higher 
order terms are indeed small. The NFE approximation has often been justified in a context far wider than linear 
response alone by appeal to the fact that for a crystalline solid, the structure dependent reciprocal lattice vectors 
are typically near the pseudo-potential zero-crossing fco with the inference that the net scattering is small [0. This 
important effect stems from the confluence an atomic and a structural length scale; the interference effect wc discuss 
is complementary, but has a different physical origin, namely an interference between intrinsic atomic and electronic 
lengthscales. Once again, the clear exception is the singular case of a point-charge (i;P^(fc) ~ 47re^/fc^), i.e. the case 
of hydrogen, which has no well-defined core- length scale fco, no oscillations in the potential and thus no interference 
effect in the higher order terms. In sharp contrast to other systems, non-linear response terms are large term by term. 
In fact, the response series may not even formally converge and great care must be taken when applying concepts 
derived from linear-response theory to hydrogen (it is not a simple material) . 

As noted, the continued accuracy of linear response is important to an interpretation now proposed for electron-ion 
structure factors Sei{k) in metallic fluids, these being defined as fc-space density-density correlation functions 
Invoking the adiabatic approximation they can always be rewritten in terms of the ion-ion structure factors as follows: 

S,i{k) - 4^5//(fc), (3) 
V Z 

which defines a new dimensionless object n(fc). Electron- ion correlations can therefore be described by convolving the 
pseudo-electron density (or pseudo-atom) n{k) with the ionic correlations. The accuracy of linear response for the 
pseudo-potential in an electron-gas implies that it should now also be an excellent approximation for a determination 
of the pseudo-atom density. For simple liquid metals Sii{k) is very well approximated by the Percus-Yevick analytic 
form for hard spheres by specifying a single parameter, the packing fraction ry, which is close to rj ^ 0.46 for most 
simple metals near melting [ p^ . Using this in (|^), we compare our approach in Fig. 2 to the full ab-initio Car- 
Parrinello calculations of de Wijs et al [|6j. The correspondence is striking, especially when we note that the 
parameters rj and Rc are a priori set by other physical properties (no fitting is necessary). 
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FIG. 2. The electron-ion structure factor Sei{k) and related electron- ion correlation function gei{r) for Mg and Bi: 
Car-Parrinello results of de Wijs et al (solid line) vs. the simple linear-response approach augmented by a hard-sphere 
approximation (dashed line). Panel (a) shows Sei{k) and panel (b) shows g^iir) for liquid Mg. Panel (c) shows Sei{k) and 
panel (d) shows Qeiif) for liquid Bi. For Mg the parameters (taken from the literature) are: rs=2.66ao, i?c=1.31ao and for Bi 
the parameters (taken from the literature) are: rs=2.25oo and _Rc=1.15ao. Both have a packing fraction ?7=0.46 (note that for 
the Qe-iif) the region inside the core radius is not physically significant). 

Besides a semi-quantitative description of electron-ion structure factors, this linear response theory now provides 
an important qualitative insight into the form of the electron-ion structure factors JTzt . The pseudo-atom density 
nik) is typically largest for small k and rapidly declines for larger fc, while the near classical ion-ion structure factor 
Sii{k) follows an inverse behavior; it is small for small k. Together with the product form (^) this implies that the 
shape of the electron- ion structure factor Sei{k) is determined primarily by the the position of the zero-crossing ko 
of n{k) with respect to the first maximum kp of Sii{k). If kg < kp, then Sii{k) selects (or filters) the negative part 
of n{k) and Sei{k) takes a form similar to that of Mg (Fig. 2 (a)). Conversely, if fco > kp, then the ion-ion structure 
factor selects (or filters) the positive part of n{k), and again, Sei{k) takes a form similar to that of Bi (Fig. 2 (^ 
(c)). Since Xi(fc) is positive definite, the zero-crossing in linear response occurs at kg. The large slope of n{k) near 
the zero-crossing then implies that non-linear corrections must have a small effect on the location of the zero-crossing, 
and together with the expected accuracy of linear response this implies that fco ^ fco- As mentioned earlier, for most 
metals, fco is just a little less than 2kF, and the latter's ratio to kp is well known: for small valence {Z < 2), 2kF < kp-, 
for large valence {Z > 3): 2kp > kp [Q. This accounts in a straightforward way for the two separate forms found 
by deWijs et al ||l^: For Mg, fco < kp {Z = 2), which belongs to the lov^r- valence class of electron- ion structure 
factors. For Bi, fco > kp [Z = 5) and we may refer to this as the high valence class of electron-ion structure factors 
p9[ . Generally ions of valence Z <2 belong to the low valence class while ions with valence Z > 3 belong to the high 
valence class. Ions with valence Z — 3 typically belong to the high valence class also, although they may characterized 
by a crossover form . The analytical approach above can easily be extended by using the modern theory of liquids 
to obtain improved ion-ion structure factors pO| , but to include second order contributions to the pseudo-atom n(fc) 
necessitates not only second order electron response, but also contributions from ion-ion triplet structure. This can 
also be carried out with concepts from the theory of classical liquids. [ pO[ 

These observations have a potentially useful experimental consequence: the principal features of electron-ion struc- 
ture factors can be measured by exploiting the differences between x-ray scattering, which probes the density fluc- 
tuations of all electrons, and neutron scattering which generally probes fluctuations of the nuclei. pl| |. X-ray mea- 
surements are usually interpreted using a free atom form factor, while our analysis suggests that for liquid metals, 
they should be interpreted with the pseudoatom as a form factor. When this is taken into account, a small difference 
between x-ray and neutron scattering determinations of the ion-ion structure factor should emerge. This difference is 
largest for metals with a high ratio of valence to core electrons. For Li (1:2) or Al (3:10), we predict a 2% difference 
at the first peak of the structure factor [E2j, but the largest effects are expected for Be which has the highest ratio 
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of valence to core electrons (1:1) and for which the difference could be as much as 7%, well within experimental 
range. In addition, Be may straddle the two classes {ko is near kp), which means that small differences in ko with 
respect to fco rnay lead to significant, qualitative differences in Sei{k), making it a particularly interesting candidate 
for illuminating nonlinear effects. In a similar way we now anticipate that higher order effects can be revealed in 
partially covalcnt liquid metals, silicon and gallium being examples. The arguments presented suggest that these 
should become relatively less important upon an increase in density (via pressure). 

The arguments and associated analyis above therefore provide a physical basis for understanding why linear response 
theories in dense electron systems generally perform so well. The accuracy of linear response is demonstrated for fluid 
metals by a simple analytical linear-response theory augmented by a hard-sphere approach to classical electron- 
ion structure factors, which already gives semi-quantitative accuracy. It suggests that there are two main classes of 
electron-ion correlation functions, one for high and one for low valence metals. Finally is it suggested that experimental 
advances in x-ray and neutron-scattering may be poised to provide measurements of these electron-ion correlation 
functions, and hence on the interactions themselves. 
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